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\Q • Abstract. We briefly recall the physical background of the transport of 

angular momentum and the mixing of chemicals inside stellar radiation zones 
and its importance for stellar evolution. Then, we describe its present model- 
ing, its successes and its weaknesses. Next, we introduce the new theoretical 
developments that allow us to treat in a self-consistent way the effect of the 
Coriolis force on the low-frequencies internal waves and its consequences for 
the transport processes. This research is aimed at improving the modeling 
of stellar interiors in the perspective of future astero and helioseismology 
missions such as COROT and GOLF-NG. 
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1 Rotational mixing 

• In standard models of stellar interiors, radiation zones which are convectively 
| stable are postulated to be without motion other than rotation. But various 
■ observational results (i.e. surface abundances, helioseismology) show that these 
04 | regions are the seat of mild mixing. The most likely cause of such mixing is 
\Q • the differential rotation, which drives a large scale meridional circulation due to 
thermal imbalance and angular momentum transport and gives probaly rise to 
t^* ■ shear instability (cf. Mathis, Palacios & Zahn 2004). This process has been 
described in a self-consistent way, namely taking into account the transport of 
angular momentum which modifies the rotation profile (Zahn 1992, Maeder & 
Zahn 1998, Mathis & Zahn 2004). Series of models have been built which include 
these processes, called rotational mixing (cf. Talon et al. 1997, Meynet & Maeder 
C$ ' 2000), and for massive stars, they agree rather well with the observations. 

However, rotational mixing alone cannot explain the almost uniform rotation of 
the radiative interior of the Sun, as revealed through helioseismology, and one 
must conclude therefore that another, more powerful process is operating, at least 
in slow rotators. The most plausible candidates are magnetic torquing (Garaud 
2002) and momentum transport by internal waves (Talon et al. 2002, Talon & 
Charbonnel 2005). 
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Concerning the waves, their treatment presents two major weaknesses. The 
first one is our crude description of their generation by turbulent convection. The 
second one is that our present description does not take into account the action of 
rotation on the waves. In this work, we have undertaken to improve the modeling 
of the transport by internal waves by introducing the effects of the Coriolis force. 
Indeed, the low-frequency internal waves which are responsible for the deposit or 
the extraction of angular momentum (cf. Talon et al. 2002) are strongly influenced 
by the rotation because the frequencies of the waves are of the same order that 
the inertial frequence, 20. Thus, internal waves become gravito-inertial waves (cf. 
Dintrans & Rieutord 1999) and we have to treat the action of the rotation on the 
waves and their feed-back on its profile. 



2 Transport of angular momentum by gravito-inertial waves 

The velocity field with respect to an inertial frame is V (r, t) — r sin QQ (r, 9) e v + 
u{r, 6, (p, t), where f2 (r, 0) and u are respectively the differential rotation and the 
velocity of the waves, r, 9, ip are the classical spherical coordinates and t is the 
time. We solve the momentum equation 



D t u+ 



2<>e ; \ (7+ ;-hin Oii ■ V<»e = - — — +— g where D t = d t + 0,d v , (2.1) 



the continuity equation D t p + V • (pu) = (2.2) 

and the energy equation which is given in the adiabatic case by 

"dlnp 1 dlnP" 




dr I\ dr 



0. (2.3) 



g is the gravity; p and P are the density and the pressure and p and P their 
respective eulerian fluctuations. We treat the case of weakly differential shellular 
rotation: 

fl (r, 6) = O (r) = H s + STi (r) where 5Q{r) <:Ti s . (2.4) 

fi s is the mean rotation rate of the region where the waves are excited by turbulent 
convection (namely the border between the radiation zone and the convective core 
or envelope which is studied), while Sil (r) is the residual differential rotation. We 
introduce the characteristic frequencies of the problem: 

er (r) = er + mQ (r) = a s + mSfl (r) where a s = <7o + mft s ; (2.5) 

a (r) is the local frequency "seen" by the Doppler shifted wave, and <7o and a s 
are the frequencies of the waves respectively in an inertial frame and in the coro- 
tating frame. We ignore the perturbation of the gravitational potential (Cowling 
approximation) and filter out the effects of the centrifugal force, which have to be 
taken into account in rapid rotators. 
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We are interested here in the regime where both 20, <C N and a -C N, 
and we make the traditional approximation, which consists in neglecting the 
horizontal component of the rotation vector in the momentum equation. We 
can then separate the variables in the partial differential equations which gov- 
ern the problem, as illustrated here for the vertical component of the velocity 
field: u r = J2k, m u r;k,m (r) &k,m {x; v) exp [irrup] exp [ia t] where v = 2Q s /a s and 
x = cos 9. The radial functions u r; k, m obey the following vertical differential 
equation 



dr 2 



+ 



AT2 



1 



A fe , m (v) 



^k,m = where ^k,r 



P 1/2 r 2 u r , k , m , (2.6) 



while the horizontal functions Qk,m (x; v), which are called the Hough functions, 
are the solutions of the Laplace equation 



_d_ 

dx 



l-x l 



1 



1 + v 2 x 2 



-mv- 



6fc,m = -Afe.m {v) 6 



k.rn i 



1 — v 2 x 2 dx J 1 — v 2 x 2 \1 — x 2 1 — v 2 x 2 

(2.7) 

which satisfy the relevant boundary conditions. Those equations have the same 
structure that in the case where Coriolis force is ignored. The difference is the re- 
placement of the classical horizontal eigenvalue of the spherical harmonics, 1(1 + 1), 



by the horizontal eigenvalue related to the Hough functions, 



the vertical wavenumber ky. k m 



>). Note that 
now depends 



N 2 /a(rf-l) -A k , m (v)/r 2 
on rotation, and therefore also the radiative damping. The second equation de- 
scribes how the horizontal functions are modified by rotation. 

Next, using the continuity equation and the quasi-adiabatic approximation (cf. 
Press 1981) to describe the damping acting on the waves, we get the following 
equation for the transport of angular momentum in the hydrodynamical case: 



4<™ 



1 8 r (pr 4 W 2 ) + ^8 r (pv v r 4 d r Q) - -^d r [r 2 Tj (r)] 



5r 2 



(2.8) 



The terms in the right-hand side of the equation are the divergence of respectively 
the mean flux of angular momentum advected by the meridional circulation, the 
mean viscous flux related to the shear-induced turbulence and the mean flux trans- 
ported by the waves, Tj, which is given by: 



1 



:N, 



47rr^Fj (r)=-r c fc-VcV?-f 



N c E \ l u{")] I 

E E ■ 

1=0 m=-l 



1(1 + 1) 



(21 + 1) 




E { i V tn M] 2 ex P Hfc,™ (r, SU (r) ; „) ] 

{k\A k:m (v)>0} { A k,m \ V ) 



where 

r k , m (r, flj (r) ; v) = A* m („) £ {^^7y - o 2 (/)] " } ^ . (2.10) 
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r c is the radial coordinate of the top of the radiation zone, p c and N c are the 
density and the Brunt- Vaisala, frequency. V c , u c and k c characterize the largest 
convective eddies which are at the origin of the generation of the waves. The I 
are the order of spherical harmonics which describe turbulent convection on the 
sphere. Here, we must underline that we have taken the simplest description of 
turbulence, namely Kolmogorov's law (cf. Garcia Lopez & Spruit 1991, Zahn 
et al. 1997). V\ m (y) is the integral related to the projection of the convective 
velocity field on the Hough functions while J\-,k,m {v) is related to the mean flux of 
angular momentum transported by waves on an isobar. Finally, K is the thermal 
diffusivity. 

Therefore, the Coriolis force is acting on the structure of horizontal eigenfunc- 
tions of the waves but also on their radiative damping. Indeed, the values of 
Afc,m (y) are greater for retrograde waves than for prograde ones (cf. Lee & Saio 
1997). Thus, the properties of extraction and deposit of momentum are modified. 
In the same way, one should also note that the values of V\ m {v) are different for 
prograde and retrograde waves. Therefore, Coriolis force modifies the respective 
generation of retrograde and prograde waves by turbulent convection. 

This treatment is the first one that allows a coherent treatment of the trans- 
port of angular momentum by gravito-inertial waves. We first treat the action of 
the Coriolis force on the waves and next their feed-back on the residual differential 
rotation. 

Work is in progress to implement the new equations in existing stellar struc- 
ture codes (STAREVOL & Geneva, CESAM) to model the evolution of rotating 
stars and to improve the modeling of generation of internal waves by turbulent 
convection (spectrum and associated flux). 
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